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Abstract 

We attempt to give a holographic description of the microscopic theory of a BCS supercon- 
ductor. Exploiting the analogy with chiral symmetry breaking in QCD we use the Sakai- 
Sugimoto model of two D8 branes in a D4 brane background with finite baryon number. 
In this case there is a new tachyonic instability which is plausibly the bulk analog of the 
Cooper pairing instability. We analyze the Yang-Mills approximation to the non-Abelian 
Dirac-Born-Infeld action. We give some exact solutions of the non-linear Yang-Mills equa- 
tions in flat space and also give a stability analysis, showing that the instability disappears 
in the presence of an electric field. The holographic picture also suggests a dependence of Tc 
on the number density which is different from the usual (weak coupling) BCS. The flat space 
solutions are then generalized to curved space numerically and also, in an approximate way, 
analytically. This configuration should then correspond to the ground state of the boundary 
superconducting (superfluid) ground state. We also give some preliminary results on Green 
functions computations in the Sakai - Sugimoto model without any chemical potential. 

Keywords: Holographic QCD, Sakai-Sugimoto Model, AdS/CFT Correspondence, 
Superconductivity 



1. Introduction 

The application of AdS/CFT techniques to strongly correlated systems in condensed 
matter situations is a very promising development [1, 2, 3, 4, 5, 6, 7, 8, 9]. This is because 
there are genuinely strong coupling fixed points in condensed matter systems in contrast to 
particle physics. Even the so called strong interactions of particle physics are described by 
an asymptotically free theory, QCD, which has a fixed point at zero coupling. Of course 
the coupling constant of QCD is large at low energies and it may be that to a large extent 
these strong coupling regions dominate the dynamics in certain phenomena. In that case it 
is plausible that one can borrow results from N=4 super Yang-Mills also at strong coupling. 
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where indeed one can apply AdS/CFT. Presumably the agreement between the small values 
of ^ seen in quark gluon plasma in heavy ion collisions and in N=4 Yang-Mills has to be 
understood on these lines [10, 11]. 

In applying AdS/CFT techniques to condensed matter physics or QCD one can take a 
phenomenological approach where one appeals very strongly to universality. This means that 
the details of the model in the UV are unimportant for low energy phenomena. Many recent 
papers have analyzed this idea and provided very useful insights in understanding earlier 
calculations [12, 13, 14, 15]. This may be qualitatively true in many situations. Thus there 
are some gravity descriptions of QCD-like theories where features such as confinement or 
chiral symmetry breaking may be seen. They can be said to be in the same universality class 
as QCD. They are not known to be reliable quantitatively (and not in the UV regime either). 
Thus the glueball spectrum is not in very good agreement with lattice results although there 
are some generic features that are the samc[16, 17]. But it may even be possible to do 
this in a quantitative way by introducing some parameters and fitting to experimental data. 
Thus in standard field theory QCD - chiral Lagrangians are motivated by this idea. This is 
also the spirit of the Landau-Ginzburg analysis in critical phenomena in condensed matter 
physics. In the context of AdS/CFT there have been some approaches to QCD with this 
philosophy [18]. 

However in condensed matter systems there are models that one takes a little more 
seriously. There are models such as the Ising model or Hubbard model, that one would hke to 
solve as a fully quantum theory. This means that while the model may be an approximation 
to the condensed matter physics, the model itself as a mathematical system is taken fully 
seriously. One cannot do this for instance in Landau-Ginzburg theory or chiral Lagrangians - 
both of which are non renormalizablc. Pcrturbative non renormalizability implies that order 
by order one has to extend the model by adding a large number of additional interactions and 
consequent additional free parameters. The model then has to be understood as an effective 
low energy theory, with a cutoff. Only then does it make sense as a quantum theory. 

This sort of distinction applies to AdS/CFT models as well. The models that start with a 
bulk gravity theory and where one does classical gravity calculations to extract physics of the 
boundary belong to the Landau-Ginzburg/chiral Lagrangian class. One cannot take them 
seriously beyond the large N limit, because the quantum gravity theory is not well defined. 
Thus one does not take them seriously as a mathematically well defined system beyond 
leading order. Furthermore the boundary theory has no independent definition. Very often 
we do not know the action or even the field content. 

On the other hand there are models where one begins with a string theory background, 
where in principle the theory is defined at the quantum level. Some examples are the models 
for QCD first studied by Witten [19] and the more recent versions of it such as the Sakai- 
Sugimoto model, where flavor D-branes are added and give fundamental fermions - quarks 
[20, 21, 22, 23, 24, 25, 26, 27, 28, 29]. Although these theories are non conformal, there 
is an underlying AdSj geometry and hence one can think of these theories as deformations 
of some conformal field theory. Apart from being well defined quantum mechanically, they 
have the advantage that the field content of the boundary theory is precisely known. One of 
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the nice features of the Sakai-Sugimoto model is that it gives a geometric picture of chiral 
symmetry breaking. Superconductivity in strongly coupled M — 2 Super Yang-Mills theory 
with fundamental matter was studied via AdS/CFT in [30] 

It is well known that the instability of Cooper pairing in superconductors is very similar 
to that of chiral symmetry breaking [31]^. This analogy motivates us to study the Sakai- 
Sugimoto model to understand strong coupling BCS superconductivity. The Sakai Sugimoto 
model has the feature that supersymmetry is broken and the lightest modes in the boundary 
are fermions in the fundamental representation, rather than scalars or some supermultiplet 
of particles (as in some of the D3-D7 models). The fact that the light particles are fermions, 
makes this realistic as one can be assured that what is condensing is a bound state of fermions 
rather than some scalars. Thus one is actually seeing the Cooper pairing. This goes beyond 
the Landau Ginzburg description of the theory as a C/(1)-Higgs system. 

BCS vs Chiral Symmetry Breaking: 

There are two important differences between the Cooper pairing instability in BCS theory 
and chiral symmetry breaking. Both these must be kept in mind while using the Sakai- 
Sugimoto model. One is the presence of a Fermi surface in the BCS situation. Chiral 
symmetry breaking as studied by Nambu-JonaLasinio (NJL) takes place in the vacuum. 
This difference shows up in the two gap equations. The NJL gap equation is a one loop 
tadpole vanishing condition. The one loop fermion contribution is balanced against a tree 
level constant and gives: 

g I d?k -^=^= ^ 1 (1.1) 
where m > is the sign of chiral symmetry breaking. The analogous BCS gap equation is 

g / d^k , ^ 1 (1.2) 

where = 

If we put m = in (1.1) and solves for g , one finds g = g^^ ■ For g > g^r one finds 
m > . This is the chiral symmetry breaking phase. In the original NJL calculation 
fi'cr ~ ^ • Contrast this with (1.2). Only states above the Fermi sea contribute and the 
criterion for the instability is A > . But note that A is measured from the Fermi surface. 
If we set A = there is an IR divergence as A; — > /c^ , which in turn is because the measure 
d^k ~ ATrkpdk does not vanish. Thus for arbitrarily small g one can find a solution with 
A > . In fact one can write (1.2) as 

dn r= , 1 
9 -r de- 



de Jo Ve^ + A2 



For g ^ « 1 , the solution is A ® . On the other hand for ^ >> 1 , we have 

^ ~ ^c9' ^ ■ The dependence of the gap A and hence also the critical temperature on the 



^Tliis is obvious from the title of the paper of Y. Nambu and G. Jona-Lasinio: 
■Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity" . 
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parameters {g, ^) is very different. In weak coupling tliere is a non perturbative dependence. 
In strong coupling (which is what the holographic calculation gives) the temperature is 
proportional to some power of the number density. 

Thus in the presence of a Fermi surface there is always a BCS instability.^ This points to 
an important modification of the Sakai-Sugimoto model. We need to include a chemical po- 
tential for "baryon" number i. e. we need a finite number density of baryons in the boundary 
theory. 

The second important difference is that it is a chiral rather than vector symmetry that 
is being broken in the NJL model. Chiral symmetry is broken as soon as the fermion has a 
bare mass, so in this sense it is a symmetry that is easily broken explicitly. In fact the Vafa- 
Witten theorem shows that vector symmetries are not broken spontaneously in QCD-like 
theories when the bare quark masses are non zero [33]. 

In QCD and also the Sakai-Sugimoto model the vector symmetry is unbroken. With 
zero bare mass for the quarks, the unbroken phase (with two flavor D8-brane pairs) has 
SU{2)l X SU{2)r = SO{A) . The quark bihnears ^qj form a four-vector under this 5'0(4). 
When one component gets an expectation value what is left unbroken is by definition the 
vector SU{2) - also called strong isospin. A basis can be chosen such that the Goldstone 
bosons are qaq - the pions. In this situation the question of whether this condensate could 
have broken the SU (2) vector is meaningless. In fact if we find that one of the components of 
the pions condense due to quantum effects, this corresponds to a realignment of the vacuum 
and one can perform an SO{A) rotation so that what is left unbroken is still an SU{2) . If 
there is a bare quark mass then there is a well defined notion of which is the vector symmetry. 
This is the situation studied by Vafa-Witten. Also the presence of electromagnetism changes 
the situation, because there is a preferred axis defined by the charge of the pions and this 
breaks the strong isospin. Thus if due to external field effects the (7r+) 7^ then we cannot 
rotate this away and we have a breaking of C/(l)em- ^ 

Thus if we are to spontaneously break a vector symmetry then the Vafa-Witten theorem 
tells us that it cannot be in the vacuum. We must have a finite number density of fermions. 
This again points to the same requirement. The Sakai Sugimoto model has to be modified: 
we need a chemical potential and a finite number density. Such modification have been 
studied, see [34, 35] for this model and [36] for other models. 

Baryon number corresponds to a U{1) charge corresponding to a gauge field on the D8 
brane. To probe this we need a field on the D8 brane that has non- zero baryon number. We 
can achieve this if we have two D8 branes, so that the U{1) can be embedded in SU{2) . 
The charged gauge fields ("w^") or charged scalar can then be used as the fields dual to 
the baryon number violating charged condensate. If the branes are separated this would 



^This can be also be seen using the Renormalization Group formalism where the BCS perturbation shows 
up as a relevant operator [32] . 

''Note that if (tt^) 7^ , when the quarks are massless, this can be rotated away. Thus if u\uR-\-d}^dR = v\_ 
(the isosinglet) and u\ur — d^dfl = V2 (the neutral component of the isotriplet), then an SU{2)l rotation 
that commutes with U{l)em can rotate one into the other. 
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correspond to a Georgi Glashow model where the SU{2) is broken to a U{l)s. (The scalar 
field on the D8 brane, corresponding to its transverse fluctuations, becomes the adjoint scalar 
fleld of the Georgi Glashow model.) The generator of this U{1) will be called deflned later 
in this section. This makes classical computations reliable and so we will implement this. 

By an AdS/GFT-type dictionary one would expect that the boundary value of the 
charged field would correspond to the value of the condensate. ^ 

In the boundary theory one has two quarks in the fundamental of N^. and a doublet of 
flavor, which we denote hj u,d . If we turn on a chemical potential corresponding to the 
3-component of SU{2) this couples to — d^L^L ■ The height of the Fermi surface will 
then be different for u and d. But in the confined phase one should define Fermi surfaces 
for color singlet fermions. The natural candidates are the baryons (for Nc odd). Thus one 
expects an excess of baryons with > ^. (If Nc is even this has to be interpreted as a 
chemical potential for scalar baryons. This is more like a BEG situation.) Then there is the 
question of what the BGS condensate consists of. Since in the bulk we are studying only the 
fields on the left branes, we are studying boundary fields involving UL,dL ■ ^ One can expect 
pairing between the ul, d^s of the form ulUl or dLdi- This would break color and give rise 
to color superconductivity [37]. Since our bulk probes do not carry color we cannot study 
this easily. Thus we focus on color singlets - mesons and baryons. One cannot make scalar 
mesons out oiuLidi alone - the scalar mesons necessarily involve both ul and ur Note 
that ud — u^dn + d\uR . The massless pseudo-scalar mesons, the pions are frozen by our 
boundary conditions on the gauge fields. The other mesons are massive. In principle they can 
still condense, but being made out of ul and ur cannot be probed by our bulk fields which 
are only on the left brane. On the other hand baryons can form Gooper pairs and condense 
as in BGS theory. Thus if is odd B = eai,a2,-,a]Ve'"L '^l^'^l "^l ' ' ^'^d can have 

(BB) 7^ 0. ^ Even though the gauge theory is strongly coupled, the coupling between the 
color singlet baryons is not so obviously strong. Thus one need not expect protons and 
neutrons to form Gooper pairs and condense in the vacuum. It certainly doesn't happen 
in real life QGD. However with a finite chemical potential and a Fermi surface for baryons 
it is possible. This description in terms of Baryons is reminiscent of the complementary 
description of color superconductivity in [37]. This is the most likely boundary dual of the 
bulk condensate that we are studying. 

From the bulk perspective adding a point source baryon charge, corresponding to the 
overall U(l), to the Sakai-Sugimoto model makes the D8 branes have a cusp geometry [34] 



^The asymptotic geometry here is not AdS. However the underlying M- theory geometry has an AdSr and 
so while the usual AdS/CFT dictionary does not apply, one may expect some modified version of this to 
hold. A similar situation involving Dl branes and an underlying AdS4 has been studied in [38]. 

^ts is defined in the next page 

''This assumes that one can treat the dynamics of the left and right brane independently. This is not the 
situation in the Sakai-Sugimoto model which is left-right symmetric. However in more general backgrounds 
with electric fields this is possible. 

^The vector mesons cannot condense without breaking Lorentz invariance. One could imagine a bilinear 
product of vector mesons condensing. 

^If Nc is even, the Baryons are scalars and can condense directly. This is like a BEC. 
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Figure 1: In the presence of a point source of baryon number there is a cusp singularity as shown on the left. 
When we have two branes and SU{2) breaking source eliarge the branes separate. In the second figure the 
mirror D8 branes are not shown in detail but are indicated by a dotted line. Their configuration is chosen 
to satisfy the force balance condition. 



(see Figure (1)). 

The coupling of this point charge to the U{1) gauge field on the brane is via a Chern- 
Simons term in the action via a point like instanton configuration. The instanton can also 
be understood as a D4 brane wrapped around the and has the same action [20] . However 
the instanton coupling requires more than one brane since it is a configuration involving 
non Abelian gauge fields. Thus we should interpret our single brane as a collection of 
coincident branes, with > 2. The only role of the extra branes is to accommodate the 
instanton configuration. Otherwise as far as the rest of the dynamics that we are interested 
in is concerned, the collection behaves as one U{1) brane. This C/(l) is the overall center of 
mass C/(l) in the U{Na). All this is exactly as in [34]. 

In our case the new thing is that we have two such D8 branes with a U{1) charge and 
Chern Simons term for each brane. Thus we have a U{2) gauge symmetry on the branes. 
When the charges on the two branes are unequal then the branes come in at different angles 
and separate. The U{2) symmetry is broken (Higgsed) to [/(I) x [/(I) when the branes are 
separated. These can be denoted as U{1)b x C/(1)3. U{1)b is the overall C/(l) of the U{2) 
and ^7(1)3 corresponds to the generator of the SU{2). 

We define the SU{2) generators (^+,^-,^3) : 

which satisfy the SU{2) algebra 

[t3,t±] = ±t±, [t+,t-]=2t3 {t+,t-} = /. 

As mentioned above the Chern Simons coupling of each brane comes from having A''^ 
coincident branes. Thus we actually have 2Nfi branes and a U{2Nd) symmetry. The U{2) 
group that we are concerned with in this paper is thus embedded in an obvious way in the 
U {2Nd) group. (The mathematical description of this embedding is given for completeness 
in Appendix B.) The brane configuration thus leaves unbroken a C/(1)b x SU{Nd) x U{l)s x 
SU{Na) group. 
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From the bulk brane dynamics also we expect a condensate, because it is known that 

when branes intersect at an angle, in general there are tachyonic excitations [39, 40, 41]. 
Thus there should be a condensate of charged fields. This should then correspond to the 
Cooper pairing instability of the boundary theory. ^° This generates a non zero electric field. 
In the presence of a non-zero electric field and charged condensate the tachyonic mode should 
get lifted and one expects a stable solution with electric field. This condensate breaks one 
more f/(l), namely f/(l)3 , in the U{1)b x SU{Nd) x f/(l)3 x SU{Nd) group left unbroken 
by the geometric description of the D8 branes. As in [34] since the SU {Nd) x SU {Na) gauge 
fields play no role in any of the dynamics we do not mention them again in the discussions 
below. 



There is a geometric interpretation of the instability that makes this clear (see Figures 
2, 3). The tachyonic mode corresponds to the tendency of the D branes to separate in 
the region around the intersection, and become parallel, because this reduces the energy. 
However, once there are open strings connecting the branes, it costs energy to separate the 
branes and make them parallel. Thus one expects stable solutions with condensates. We also 
expect to find analytic solutions in fiat space describing charged condensates. The absence 
of tachyons in the presence of an electric field is confirmed by a perturbative analysis of the 
fluctuation spectra. 

In fact there are solutions to the Yang-Mills equations that describe precisely such a 
situation: as soon as we have a non zero t-^ charge, there are solutions where fields 
condense. If the condensate extends to the boundary, we expect superconductivity in the 
boundary. We also give a proposal for calculating Greens functions in this situation. The 



'We are assuming that there are no other instabiUties in the boundary. 





Figure 2: Tachyonic instability 



7 



Figure 3: Electric field causes open strings to be formed that stabilize the configuration. 



complication here is that the brancs do not extend all the way to the black hole. Hence 
the usual ingoing boundary conditions at the horizon cannot be applied because there is 
no horizon crossing. We apply this technique to the Sakai-Sugimoto model where Greens 
functions have not been calculated. We expect that the boundary fermions being massive the 
Green's function (say, of the currents) should show a gap in the imaginary part. We find that 
the conductivity calculated using this prescription seems to have the expected properties. 
However we leave a detailed analysis to the future. Also applying to the BCS situation is 
more complicated and we have not attempted it in this paper. 

We can summarize the above discussion as follows: In the bulk we have two (sets of 
coincident) D8 branes in the configuration shown in Figure 1. There are two C/(l) charges 
di, d'i corresponding to each brane. This can be represented as 

Ml \ , , , , / 1 \ 4 / 1 \ 

( dJ=^°I + ^3t3 = cio(o , ) . (1.4) 

do is the U{1)b baryon number charge and d^ is the U{l)z (corresponding to t^) charge. 
When (is > the branes come in at different angles at the cusp Uc- For u» Uc the branes 
arc parallel but separated by a finite distance. The SU{2) is then broken to U{l)s as in the 
Gcorgi-Glashow model as soon as d^ > 0. As long as no off diagonal SU{2) fields are turned 
on the DBI action is exactly that of two independent branes and is easily analyzed. The free 
energy of this configuration can be studied for arbitrary values of dQ,d^. 

We have also seen above that in this situation, a flat space analysis shows that there is a 
tachyonic instability which wc plausibly identify with the BCS instability - assuming there 
are no other new unexpected boundary instabilities. 

To quantitatively analyze this condensate in the D4 brane background we restrict our- 
selves to dz « do and all SU{2) fields smaller than the other terms in the DBI action. In 
this case it is consistent to expand the DBI action and keep the Yang-Mills part. A flat 
space analysis of Yang-Mills also tells us that there is a stable ground state where some 
charged fields condense. This encourages us to explore solutions of the Yang- Mills equations 
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of motion in the D4 brane background. The curved space equations are then solved exactly, 
numerically and in some approximation analytically. There are solutions where the fields are 
small, validating the truncation of the DBl action to the leading Yang-Mills part. Also the 
existence of a charged condensate breaking the ^^(1)3 symmetry implies that this approxima- 
tion captures qualitatively the phenomenon that we are trying to describe. The smallness of 
the fields assures us that higher order terms can only change things quantitatively by small 
amounts - not qualitatively. Thus we have a self consistent scheme that gives us a ground 
state with a charged condensate. The f/(l)3 symmetry is broken. Thus we expect the same 
symmetry breaking in the boundary also. This is then identified with the BCS phenomenon. 

One also expects that at finite temperature the positive {massY due to temperature 
effects will overwhelm the tachyonic tendency of the intersecting brane. The negative {massY 
is proportional to the angle between the brancs, which in turn is fixed by the number density. 
So at a (dimensionless) roughly equal to the number density we expect superconductivity 
to disappear. This dependence on the power of number density is reminiscent of the strong 
coupling expression given earlier (for the gap). However we need a detailed calculation before 
anything concrete can be said. In any case it is very different from the weak coupling BCS 
expression. 

This paper is organized as follows: In Section 2 we describe the brane solution in the 
presence of charges. In Section 3 we describe the Yang- Mills solution in fiat space with 
condensate relevant for the D4-D8 case and analyze its stabihty In Section 4 we give the 
solutions (exact numerical and approximate analytic) again in the D4 background metric. 
This establishes the existence of a charged condensate. In Section 5 we study analytically 
the phase structure at zero temperature, function of chemical potential but without 

any condensate. With a condensate we give some numerical results. (The finite temperature 
analysis will be done elsewhere.) In Section 6 we give a sample calculation of a Greens 
function in the Sakai-Sugimoto model using a new prescription. We conclude in Section 
7 with a discussion of our results and some open questions. The Appendices contain a 
collection of some useful results that are helpful in obtaining the equations. 



2. Background Brane Configuration 

The ten dimensional D4, brane background metric is 



3 , , 3 



ds' = " {r,.,dy-dyP + fdO') + (^^) ' (^ + U'd^l) 



Here y°'{a — 0, 1, 2, 3) and 9 are the D4-brane directions, dill is the line element of S^, €4 is 
the volume form in 4 dimensions, and Q,^ — ^ is the volume of a unit 4-sphere. C3 is the 
dual gauge field of the D4-brane and has components along 5"^ . = ngsNJl sets the scale 
of the space time curvature. Ukk is related to the supersymmetry breaking scale, which is 

the period of the 9 coordinate, ^xl: = 27ri?r4 . 
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One can define new dimensionless coordinates: u — ^, x"" — ^ and t — ^, as well as 



where f{u) = 1 — with ukk = ^^^^ , and e*^ — Qs . 

Our configuration consists of two flavor D8 branes with a background gauge field and some 
delta function sources. (As explained in the Introduction, these two D8 branes are actually 
two sets of Nd number of coincident D8 branes.) The delta function sources correspond to 
baryons that are uniformly distributed in the {x^,x'^,x^) directions and at fixed u — Uc ■ 
Baryon number can come from D4 branes wrapped around and immersed in the D8 brane. 
They can equally well be thought of as an instanton background gauge field configuration 
of the non Abelian gauge field on the D8 brane. As shown by Sakai and Sugimoto, either 
picture yields the same value for the action. 

The Chern- Simons action is (we have integrated by parts the action C3AFAFAF and 
separated a C/(l) part for ^0) 



with J^4 dC^ = 2nNc and F A F = 87r^A^4 J duS{u — Uc) where we have assumed the 
instantons are localized at Uc- 

As mentioned in the introduction even when we have one D8 brane we have to think of 
it as a set of coincident D8 branes with the associated non Abelian gauge fields providing 
the localized instanton configuration [34, 20]. These non Abelian fields are zero everywhere 
else and play no role in the rest of the dynamics. Thus we assume following Appendix B 
that A = A^A"" + /1A°, where A" are the Gell-Mann A-matrices for SU{3) and A° is the 3 x 3 
identity matrix. The SU (3) part of the gauge field is assumed to correspond to a localized 
instanton, so that F A F — Sn'^N^ J duS{u — Uc). Aq — AqX'^ above is then assumed to 
be the f/(l) part. 

When we have two D8 branes we assume two decoupled sets of the above, each with 
a U (3) gauge group. The embedding of this U{3) x U (3) in U (6) is given in Appendix B. 
As mentioned in the introduction, these SU (3) fields play a role only in the Chern Simons 
action. They can be set to zero elsewhere. The effective DBI action is just a C/(2) non 
Abelian DBI action. 

When the branes are separated, and the non- Abelian fields are not excited, this reduces to 
two decoupled Abelian DBI actions. In field theory language, we have a.U{l) x SU{2) group 
broken to a [/(I) x ^7(1). The Higgsing of the SU{2) is by the adjoint field corresponding 
to the separation of the branes. This is just the Georgi-Glashow model. C/(l) x C/(l) can be 
denoted U{1)b x [7(1)3 as in (1-4). Thus A = AbI + A^ts. 
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We can factor out a 3- volume V3 of the three dimensional space of the boundary theory, 
by writing N4 — . ^4 is also dimensionless. In terms of these we get 

Scs = N.N^ J dtEAoiuc) = UBfaoiuc) (2.2) 
where, integrating over the Euclidean time ts with circumference /3, we have 

The instanton number A^4 is thus the number of DA branes and equivalently the number of 
baryons. The charge of a baryon is taken to be Nc- 
The action of iV4 DA branes wrapped on 5*^ is 

where we have integrated over the Euclidean time Ie with circumference j3. From the Uc 
dependence we see that there is a downward force. 

Using A^4 = 714;^, Qsll — and ^4 = ^ we get for the D4-brane action: 

{2'KYRni " z ^ ' 

2.1. Solution of DBI action with source 

The first step is to find the background D8 brane configuration. To this end we solve 
the equations of motion of the D8 brane Dirac-Born-Infeld action with a background gauge 
field on the D8 brane worldvolume. 



'>DBI 



-fjda^ [u" + u'fr" - /(9,ao)']^ - nsfaoiuc) (2.4) 



where u' = d^u , t' = daT , ub was defined above and gives the magnitude of the charge to 
which oo couples, and T — ■ Let 



Then the equation of motion for oq is 

da Ua = Ub S{u - Uc) . (2.5) 

Choose the coordinate a along the brane to be equal to u, the bulk radial coordinate. 
Then u' — 1 and r' = 9„t . 
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Case 1: ub — 



We then find in the absence of sources that Ha = do, a constant. Thus, 



2 



l + «VV2-/(a„ao)2 



dl (2.6) 



Solving for duttg we find 



^^"""^ - / {u^ + dl) ^^-^^ 

One can plug this back into the action and evaluate the Legendre transform of the Lagrangian 
density: C = C — T Ua d^ao : 



C{do) = f ^u^ + dl + u^fr'' (2.8) 

We can repeat this for t' — duT and write the Legendre transformed Lagrangian density 
jC{co,do): 

£(co,do) = -r^ (2.9) 

where D = u^f + d^u^f — Cq . When D vanishes (for some u = uq), r' = oo. This is the 
lowermost point of the D8— branc and here the vertical force due to the brane tension vanish 
because the brane is horizontal. The horizontal forces cancel if we add a mirror D8— brane 
so that we have a symmetric U-shaped configuration. In the Sakai Sugimoto interpretation 
the second set of branes provide fermions of the opposite chirality, and are referred to as 
L>8— branes. The joining of these two sets of £)8— branes represents the breaking of chiral 
symmetry. 

Naively the gauge field strength diverges at Uq. However this is a "coordinate singularity" 
- the choice a = u is not good at u = Uq because a keeps increasing monotonically as we 
proceed along the bottom of the U-shape, whereas u starts to decrease again. See Appendix 
A. In this appendix a more general analysis is done keeping the coordinate a throughout. 
Thus {dutto) = -^^^ and (St-Oq) = -^^^^ can be calculated for a given solution. The former 
is singular at uq because u^j vanishes at -u = Uo, whereas the latter is finite. Calculating the 
latter can also be accomplished equivalently very simply by choosing a = t and investigating 
the same solution. One finds: 

(a.u.r - 1^ . (2.11) 



^^We use the notation f„ = ^ 
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Note that uq is a solution of D{uo) — . So ul + d^— j-rz^ — y • This gives 



'2 



Cq 

which is finite. 

In the absence of charges, continuity of flux would require that the electric fleld continue 
in the same direction along the other brane and reemerge on the boundary. Thus dfjGo is 
continuous. On the other hand 9^ao changes sign because a — — r on the mirror brane. 
Since the equations above fix only the magnitude of 9,-ao, this is also a valid solution. The 
conclusion is that it is possible to have do ^ even when ub — 0- 

Case 2: ub 

In this case there is a jump in the value of ![„ at m = Uc- We will choose a solution where 
Ua — for u < Uc and equal to for m > Uc- 

duT = , u>Uc (2.12) 

U2j D 

where D = u^f + n\u^ f — Cq . Uc is determined by minimizing the action with respect to 
variations of Wg, subject to the constraint that 



/ 

J u, 



OO T 

duT'^- (2.13) 



L is the distance between the D8 brane and the anti-D8 brane (assumed symmetric) at 
u — OO. Note that the above constraint implies that 

dr' 

du-^- t'{Uc) = (2.14) 

(Jj tvQ 



We work with S{d) — —T du^/u^ + y j + ^^/^'^ for the D8 brane. The variation of 
Uc gives the upward force due to the D8 brane tension. In order to get the downward force 
we (following [34, 35]) use the action of a D4 brane wrapped on . The answer evaluated 
above isj^uc . This gives a force of We do not need a horizontal force balance as the 
mirror DS— brane automatically provides the correct balancing force. 

One can do the variation keeping either L fixed or Cq fixed. The latter is easier, since we 
just work with the Legendre transformed action: 



du 
u^f 



u^f + nlu^f-cl (2.15) 



Varying Uc just sets the integrand C{cQ,nB) equal to the force due to the D4-brane. 
One can solve for Uc in terms of cq^ub- The result for cq = IQQ^ukk = 1 is plotted in the 
figure. 
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If one substitutes for cq in terms of r' one can manipulate the D8-brane variation into 
the form: 

1 3 rz — TV 1 



T 



—3 u 



B 



2^/2 



tension 



sin a 



In this form it is easy to see the physical interpretation. The factor {sin a) gives the vertical 



component of the force. One also finds siv?ac 
to the horizontal aX u = Uc ■ 



ul+ng 



^ where ac is the angle with respect 



50 ino 150 

Figure 4: Uc versus do{— ub) with cq — 100 
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2.2. SU{2) Breaking 

Our aim is to introduce a finite "baryon number" corresponding to the U{1) subgroup 
of SU{2) generated by ^3 . When this is done there are configurations of the type shown 
in the second diagram in Figure 1 where the two branes separate. The branes are parallel 
at large u, but separated by a finite distance. This Higgses the SU{2) to U{1) as in the 
Georgi-Glashow model. 

The question is whether there will be a condensate that spontaneously breaks this U{1) 
as happens in the BCS case. Note that on the boundary this is a global symmetry, so in 
this sense it is like superfluidity. However what is non trivial is that a fermion bilinear must 
condense in order for this symmetry to break. This is more like BCS than condensation of 
a bosonic field. 

Uc has a (rather complicated) dependence on do{= ub) and cq . Also, the angle of the 
brane depends on n^. We have to choose values so that Uc has the same value for both 
branes. Thus to begin with ci = C2; di = d2 . We then turn on di — d2 , but change Ci, C2 
such that Uc remains the same. This can be done numerically. 

Thus if we turn on a charge corresponding to ^3, which is the difference n(3) = n^i — 71,52 , 
then we expect that the two branes will subtend different angles at Uc ■ Thus the situation is 
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as shown in Figure (1). For small values of ub and the angle etc , Oici^di ~0ici,d2 ^ ''^bi —^32 — 
n(3) . 

When two branes meet at an angle, in general there is a tachyonic mode. This mode 
corresponds to the tendency of the branes to smoothen out the angle and straighten out - 
see Figure (2) . In the presence of the electric field the situation is somewhat different. The 
presence of the electric field implies that there must be open strings which connect the two 
branes and carry the electric flux. The tension of these open strings then in principle can 
prevent the separation of the branes. In fact a perturbative stability analysis shows that 
there is no tachyon. 

In general one expects the relation between the slope ccc and the electric field to be 
modified by the condensate. The condensation corresponds to 0=*= and in general also w'^. 

The background configuration breaks SU(2) . The Non Abelian DBI action [42] is well 
defined as long as the background is diagonal, which means it lies in the f/(l) x [/(I) subgroup 
of U{2). Otherwise there are ambiguities in the definition of the determinant . There is the 
symmetrized trace prescription that is often used, though in general it is known that there 
are modifications at the 0(q;')^[43, 44]. 

When the SU{2) fields are small one can study this problem in the Yang-Mills approx- 
imation to the Dirac-Born Infeld. In Section 3 we give an exact solution of the Yang-Mills 
equations in fiat space. This solution is not the most general. It has 0"*" = , although 
w'^ . As a consequence the asymptotic slope is still fixed by the electric field. 

In the presence of the D4 D8 background and the supersymmetry breaking metric function 
f{u) , we find an approximate solution of the Yang-Mills equations by neglecting some terms 
in the equation of motion. We find that at least when Uc » RAdS the neglected terms are 
small and we expect our solution to be a good approximation. This is described in Section 
4. To do better than this one needs a numerical solution. We give some exact numerical 
solutions in Section 4. The solutions where the SU{2) fields are small a posteriori justify the 
Yang-Mills approximation and are thus very close to the exact solutions of the DBI action. 

The conclusion is thus that as soon as we turn on a charge (corresponding to t-^) we have 
a solution where there is a charge violating condensate. This must be a reflection of the fact 
that from the boundary theory viewpoint, as soon as we have a Fermi surface (i.e. a non 
zero chemical potential) there is an instability towards a ECS phase. 

As a caveat one should note that this analysis does not take into account the tachyon 
condensation that represents the fermion mass generation. However this should not affect 
the BCS instability, which only requires gapless excitation above the Fermi surface. The 
existence of this is not affected by the fermions being massive. 

In the next section we analyze the stabihty of this solution in flat space using Yang-Mills 
as an approximation to Dirac-Born- Infeld. 



^^As mentioned earlier this conclusion is contingent on the absence of other instabilities in the boundary 
that give rise to the same kind of symmetry breaking. 
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3. The Yang-Mills Approximation 

We will assume that the fields on the DS brane depend only on {t,u) , thus keeping 

translational invariance in (x^,x^,x^) . The S'^ coordinate dependence is also not assumed. 
The D8 branc has only one transverse direction, r, which we call here. The index fi on the 
Yang-Mills gauge field has nine values. However, we set = along the four directions 
of the S'^ . For the stability analysis we concentrate on Aq, A^ and (p. For other purposes, 
such as calculating correlators of the boundary theory, we also use A^i . 
We define the SU{2) generators {t+,t-,ts) : 

-^{11) (3.1) 

which satisfy the SU{2) algebra 

[t3,t±] = ±t±, [t+,t_]=2t3 {t+,t_} = I . 

If we write A = A+t+ + A't. + A% then Tr{AA) = l[{A^f + AA+A'] . 

In this basis we define the Yang-Mills field ^4^ = ^^^+ + + ^^^3 and, similarly, an 
adjoint scalar field = + (p~t_ + 0'^t3 where 0^ is real and 0^ are complex conjugates 
of each other. The covariant derivative is defined by 

D^(l> = d^(l) + i[A^, 0] 

and the field strength by 

F^i/ = d^Ay — dyA^ + i{A^, AJ\ . 
We can write the Yang-Mills and scalar field action as 

S ^ -^-Tr[F,,F^^]-Tr[D,(t>D^(t>] . (3.2) 

3.1. Solutions of the Equations of Motion 

We first study the case where the background space time is fiat. Later, we will study the 
case where the background space time is curved. 

3.1.1. Solution with = 

The simplest non trivial solution is of course the Higgs phase where (j)^ — const. The 
resulting unbroken symmetry is just U{1) . In the brane language this corresponds to a 
constant separation of the DS brancs in the r direction. 

The configuration of intersecting DS branes corresponds, in the Yang-Mills description 
to having = qut^ where q is the slope. Clearly this satisfies the equations of motion 

where D^i^*) = 9^(*) + ilA^, (*)] . For the given configuration Du4> is constant and A^ is 
zero, so it clearly satisfies the equations of motion. 

When we embed D8 branes in DA brane background, as the discussion in the previous 
section shows, there are solutions that interpolate between these two solutions. However this 
requires a non zero ^o- 
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3.1.2. Solutions with A 7^ 

We can have a configuration where = —Eu , in addition to the above . This 
corresponds to a constant electric field E and also clearly satisfies the equations of motion. 
It is this solution (near u = u^) that interpolates with the = const solution near m = 00 , 
in the DA background. 

Solution I: 

We choose a gauge where = = . By imposing this everywhere including the 
boundary, we are freezing the Goldstone modes (pions) [23]. We consider time independent 
configurations and also set the charged scaler to zero. We list below the various equations, 
after the choice of gauge : 
The equations are: 
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dlAl - Aw+w-Al = (3.3) 

9^0=^ - Aw+w-(t)^ = (3.4) 

d^{w+Al)+w+duAl = Q (3.5) 

{{Alf - {<t>^f) w+ = (3.6) 

duiw^cj)s) + w+duh = (3.7) 

Eqn (3.5) implies that {Aq)'^w^ = c, where is c is a constant. Clearly 0^ = Aq is required 
too. Finally plugging into (3.3), Gauss law, we can integrate to get 



^li-) = = V + 7; w- = w-^ (3.8) 

The D8 brane has the expected profile - the intersection region is smoothed out, but it does 
not straighten out fully, and there is a condensate of representing open strings connecting 
the two D8 branes. This solution was also studied in [46]. 
Solution II: 

A solution can also be obtained if there are two adjoint scalars 0i and 02 , with the 
action given by 

S = -^Tr[F^,Fn - Tr[D^(PiD^4>{\ - Tr[D^<i>2D^<i>2] + Tr[0i, 4>2f • (3.9) 



We define 



= (0? + i(/.^)t3 + (0+ + i(j)^)t+ + (0r + i^2)t- 



^^Note that gauge invariance of the action implies the following relation between the equations: Dii-§j^ + 

ie[0, ^] = . This means that after the gauge choice, although we seem to have an over determined system 
of equations, the equations are not all independent. 
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and, in an obvious way, 



(x')*t3 + (X+)*t+ + (X-)*t-=01-^02 




■+ 



Of course if we set x = we have the same solution as before. Furthermore, if we set 
one of the scalar fields, such as 0i = , then the system is again the earlier one and we only 
have an analytic solution where 02 is also zero. So we try to set a different set to zero. One 
can try for instance to set 0f = = 0^ 

The equations and details of the solution arc given in the Appendix F. 

When all the dust settles we find that can be nonzero: either a constant or linear in 
u. But this changes fairly dramatically the behavior of - it becomes exponential rather 
than linear. If we take X'^ = Xo = constant then Aq can be solved for in closed form: One 
finds using the same methods as earlier 



Note that duAl{0) — , and thus the electric field vanishes at w = as required by symmetry. 
3.2. Stability 

In this section we study the stability of the solutions given in the last section. The stability 
of the solution with linear 0^ has been analyzed by [41, 45, 46, 47] and as mentioned earher 
a tachyon was found. This tachyonic instabihty is illustrated in Figure 2. 

The instability would continue till the initially intersecting brancs straighten out com- 
pletely, if we take this linearized picture seriously. However what actually happens is more 
complex. The condensing tachyon is a charged field. Once there are charges one has to take 
into account electric fields. The situation when there is an electric field present is described 
in Figure 3. Since the electric fields point in opposite directions on the two branes (we are 
considering the U{1) generated by ^3 of SU{2)), as the brancs separate open strings need to 
be formed to carry the electric flux. These open strings between the separating branes have 
tension, oppose stretching, and do not want to stretch beyond a point. This stabilizes the 
configuration. 

This intuitive picture is consistent with the linearized stability analysis given below: One 

finds that {massy — ~f (2n + l) when q'^ > E'^ > , and {massy — (2n — 1) l^l when 

£■ = , with n = 0, 1, 2, • • • . Thus, the lowest {massy oc {q^ - E'^)^ > when 7^ , and 
the lowest {massy < when E = . 

The Nf = 2 number of D8 brancs in our case are in a background space curved by the 
D4 brancs. Nevertheless near the cusp at u = Uc , where these two DS brancs intersect, the 
local physics can be studied using the flat space model. This is what is done below. 




{u) = ^/ax^ cosh (xo u) 



(3.10) 
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The equations for small fluctuations is derived in the Appendix G. 
Let us assume (j)'^{t,u) ~ (j){u) e"*™* and w^{t,u) — iw{u) e~*"^* . Using 0^ = qu and 
Aq = —Eu , we then get 

dl(f) + q{2 + udu)w + {m + Eufcf) = (3.11) 

{q\'^ - {m + Euf)w = q{l - udu)(p (3.12) 

Define 

b ^q -E , C^-g^ , U^U--g^ 

and the functions P{u) and Q{u) by 

P{u) = , ^^7' , . , Q{u) = qW - (m + Euf = h'u' - c . 
[m + Eu) Q[u) 

Then using (3.12) the equation (3.11) becomes 

9> - - (Q - P) = . (3.13) 

For large u , we have uP ~ , P ~ , and Q ~ . So, asymptotically for large u , this 
is a Schrodinger equation for harmonic oscillator. Let (j) = e~^'^ F where b = ^J (f' — E'^ > 
and u has been defined earlier. We then have 

dlF - {2bu + uP) 9„F + (c - & + (1 + hu,u)P) F = . (3.14) 

Letting F{u) = , where n > is an integer, and noting that P ~ for large u , we 
get from the coefficient of rt" the eigenvalue condition: c = (2n + 1) 6 , which shows that 

777,^ _ {q -E )^ _|_ 1) is always non-negative. When q^ — E^ , the spectrum is gapless. 

This is to be contrasted with the situation where E — . Observe that P ~ for 

large u if E , irrespective of how small is. But, if £■ = exactly then P ~ ^ for 
large u . This results in the (buuPF)— term contributing to the eigenvalue condition which 
is now: c = (2n + 1) b — which gives = (2n — 1) |g| since b = \q\ and c ~ n? when 
E = ^ . This condition reveals the presence of a tachyonic mode when P = . 

The new solution with ^ 0, described in the last section has the same asymptotics as 
the above solution, hence we do not expect a tachyon. Nevertheless, since q — E, there are 
massless modes and so we expect that there are continuous deformations. This is certainly 
true since the solutions are characterized by integration constants (c, d) which are free 
parameters and can be changed continuously. 

Finally, the fact that the negative (mass)^ is proportional to the angle, which in turn, 
is proportional to the number density, suggests that the critical temperature at which su- 
perconductivity is lost is proportional to a power of the number density. This is different 
from the usual BCS relation and is closer to the strong coupling ECS result given in the 
Introduction. 
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4. Solutions in Curved Space Time 

In this section we give the curved space counterparts of the fiat space equations given in 
Section 3.1. We also attempt to find a solution that is the counterpart of the one presented 
there. Our aim will be to establish in a qualitative way that the fiat space solutions presented 
in Section 3.1 have a generalization to the D4D8 background. A quantitative (numerical) 
study of these equations is not attempted in this paper. 

4.I. Equations of Motion 

Our starting point is the Dirac-Born-Infeld action, with a — u , 

Sdbi ^-TJ d'x du-^ [1 + u^fr'^ - /(9„ao)']5 - nBTao{uc) (4.1) 

For the moment we ignore the source term. We have 

[1 + u'fr" - fid^aof]'^ ^ 1 + lu'fr" - ^/(^.ao)^ . (4.2) 

The last two terms are nothing but ^g'^'^grridu'TY + |fl'°°5'""/(>u which are the world volume 
scalar kinetic term Qt-^ df^rd'^r and the Maxwell term f^uf^" written in curved space- 
time. The factor g^-T is a component of the background space-time metric in the transverse 
direction. 

For the case of two D8 branes, we need the non-Abelian Dirac-Born-lnfeld action [42]. 
However this is not very well defined although there are prescriptions that are known to be 
consistent up to some order [43, 44]. We will bypass this complication by expanding the 
Dirac-Born-lnfeld action as was done above and using the Yang-Mills approximation for the 
non Abelian gauge field and keeping the Dirac-Born-lnfeld structure for the U{1) gauge 
field . In order to get the correct normahzation of the the commutator term [Ao,t]^ we 
will start with a C/(l) theory that has time derivatives of r , covariantize and then set the 
time derivative to zero. Sec Appendix C and Appendix D also. 

Thus our starting point is: 

[{goo + gTrf^){guu + grrr'^) + fou]^ 

= [googuuil + g''grrf'){i + ^^V^r'')) + flJ 
= [googuu{{i + g''grrf'){i + ^"V^r'^) + ^°^"Vol)]^ 

= [gOOguu{{l + ^"VrT'' + ^°°5"Vol) + (1 + g''''grrr'')9''9rri')]^ 

Covariantize now r and set the time derivative to zero : f Dqt i[AQ,T] . This gives: 
Sdbi = —T / d'^x du —j= \fx — n'^^^^ao{uc) — n^^*"^"* ■ Ao{uc) (4.3) 
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where X is given, denoting the SU{2) part of r by , by 



+ (1 + ^"VrT'')rtrr 2Tr(i [Aq, 0])^ 

Note that the factor of 2 for SU{2) traces is because the ^3 generator is defined to be |t3 
where is the Pauli matrix. 

For the SU{2) part, the covariant derivative D^* = 9^ * has been introduced 

in place of the ordinary derivative, and F^^ = d^A^ — d^A^ + A^] is the non-Abehan 
field strength. Had the SU{2) fields been entirely in the diagonal ^3 direction, this action 
would have been exact. This was discussed in Section 1. To the extant that there are off 
diagonal terms, this action is not correct. It can also be shown that if the off diagonal terms 
are entirely in the antisymmetric Fqu part, then the symmetrized trace prescription gives 
exactly this action. (A proof is given in the Appendix C). Since the symmetrized trace 
prescription is itself known to be correct only up to 0{{a'Y) we will not belabor this point 
here. However we assume for the moment that the off diagonal terms are small and work 
with this action and expand the SU{2) part in the square root keeping the background U{1) 
inside the square root. 

Defining A = 1 + grrg^'^T"^ + 5'°°5'"''(9uao)^ , we have 

I- 5 

Sdbi ^-f J d^xdu^^/A[l + — [g^^g,, 2Tr{D^<pf 
+ (1 + 2Tr{z[Ao, <P]f + 2Tr{Fouf] ] 



+ (1 + g--gV)g''grrTr{t[Ao, 0])^ + g''g--Tr{F,^f ] ] . (4.4) 

The pre factor multiplying the Yang-Mills action has been denoted symbolically as \J —g. 

Note that if we had g'^'^ instead of g^^^ then in the limit / = 1 the scalar could have been 
thought of as another vector and the action would be the same as the fiat space action up 
to an overall factor. It is easy to see that in this case the fiat space solution given in the last 
section can be generalized very directly to curved space. 

We are interested in static and translation invariant (in the (x^, x^, x^) directions) solu- 
tions of the equations of motion. The fields thus depend only on u. The equations are the 
following. These equations assume that r' = , i.e. neglect the bending of the D8 branes. 

-du (v^^%""9„Ag) + 4 ^gg"^ {g^^w^w' + ^..0+0-) Al 



5AI 



2l ^^"Vr (0~i^n0+ - 0+/^n0-) 



bA\ 
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dS 



5S 



5S 



= 2i i^g^'g^'^w+Al) + 2i ^g'^g^^^w^d^Al - 2v^^™^,,A^0=^0+ 



5S_ 



29„ (v^^"Vri^«0+) - 2i V^g^^grrw+D^<f>' - 2y/^g'''grr{AlW (4-5) 



Here, we note the following. Let 0^ 
then follows from ^ and equations above that Q 

Kir ^TT 



It 



where 

and = w e^*^ where (p and w are real, 
? + I = constant . 

As mentioned above if we replace g^-j- by g''^'^ and set f{u) = 1 (which is a good ap- 
proximation if rt >> ukk), then one can set 0"*" = as in flat space and recover the 
same solution. Thus one would get {Aq^w'^ — \f—g (0^)^u'+ = const, which except 

3 

for the factor y/—g is exactly as in flat space. But grr = U2, so clearly one cannot do 
this. However this suggests a change of variables: Introduce auxiliary fields ib"- satisfying: 
grrDu(j)"' = g^'^D^ip"' . In all the equations where 0" enter linearly, this change gets rid 
of 0" and those equations look exactly hke the flat space case. However where 0" enter 
quadratically, one ends up with both 0" and ip"- . This happens in the gauge fleld variations. 



However, since -0^0" = (f?"^^) D^ip 



u 



Duip"" , one expects that -D„0" are very small 



compared to Duip"" for large u . It is plausible that terms involving 0" are much smaller than 
the other terms, at least for very large u . Thus our strategy will be to set these terms to 
zero as the zeroth approximation. The equations of motion can now be solved analytically 
for . In terms of ip"" one can solve for 0" . In the case where we set A = 1 , this can in fact 
be done analytically, though only a numerical solution is possible in general. One can now 
verify whether there is some range of parameters where the neglected terms are really small, 
so that we have a self consistent procedure. We do find for very large u ~ 1000 and specific 
values of the parameters, the correction terms are a few percent of the terms we keep. 

We give the substituted equations below. The terms involving 0" which are supposedly 
small are in boldface. 



5AI- ^" 



6S 



SS 



54) 



-2i 



-g {w duip'^ — w'^duip ) — 4 \/—g w'^w 



(4.6) 
(4.7) 

(4.8) 
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-- 2i du i^/^w+Al) + 2ty^ w+duAl - 2 Al(P^(P+ (4.9) 



SA-, 



ss 



+4^/^(f)+('^j}+w- - '0~M?+) (4.10) 



54)- 

-4 {w-il^^ - w+i/j-) w+ - 2y/^(Al)^(t)+ (4.11) 

Equation (4.7) can be satisfied identically by the reafity conditions: (f)"^ — ± i (t){u) , 
t/''^ = ± z V^(m) and = w{u) . The following configurations solve the equations (4.6, 4.8, 
4.9, 4.11), when the terms in bold are neglected: 



,3 DAG{uy 2Cl ^,,3^2 / X ^ 



■2 



D-iG{'aY 2C\ 



= \l ^^^^W^ + ^ , V^ii^'f w{u) = C, 



, V^9m'w{u)^C^ (4.12) 
where G{u) — G{uc) = du' ^ ^ ^ ,^ ■ The above equations imply that 



Ca C3 C4 



(4.13) 



Given tjj"- one can solve for (p"" . There are some constants of integration that can be fixed 
by choosing boundary conditions for 0'^, 0"*" . 

If equation (4.10) is to be satisfied by a finite A^ then 0^ cannot be neghgibly small. 
Thus we choose the boundary condition on 0^ so that (4.10) is satisfied with finite A^ . This 
implies that 0^ cannot be very small. It is also found that 0^-0'^ has approximately the same 
functional form as {A^y only for large u . So we restrict ourselves to large u > 100 . We also 
choose therefore to set C+ — { =^ ■^^ = 0) so that terms involving 0^-0+ are neghgible. 

Finally the various constants of integration have to be adjusted so that the terms in 
bold are smaller than the others in each equation. For example, with the boundary being 
at M = 1000 , Ca = 10"^; C3 = lO"^; Da = 10"^ C+ = 0; 0^(1000) ~ 4.47 * 
10~^; 0"'"(1OOO) = . We also set A = 1 as an approximation. 

Thus for these values, as an example, (4.10) evaluated a,t u — 900 is ~ —3.97 * 10~^^, 
and a typical large term in the equation {y/—g{AQy'w'^)\u=9oo — 5.81 * 10~^ . The same 
quantities evaluated a,t u = 100 are ~ —3.02 * 10~^^ and 1.0 * 10~^ respectively. This shows 
that the equations are satisfied to a very good accuracy. 
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As mentioned in the beginning of this section, all the restrictions have to do with justifying 
perturbation theory around an analytic solution. In order to get a more general solution one 
has to resort to numerical methods. 

Some exact numerical solutions to the equations (excluding (4.10)) are given below along 
with the analytic approximations. 

Figure 9 and 10 give a comparison of some of the analytical and numerical solutions. As 
expected the agreement is very good for large u. 



Allu] 




Figure 5: Plot of A^- numerical. The second plot expands the small u region. 
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Figure 6: Plot of the charged scalar field - numerical. The second plot expands the small u region. 



A technical comment about the numerical solutions is that equation (4.10) is not auto- 
matically implemented. However gauge invariance of the action implies the following relation 
between the equations: D^^j- + ie[0, = . Thus, in this context it means that 




(4.14) 



This means that if the other equations are satisfied (4.10) must evaluate to a constant. Thus 
for some choice of boundary conditions this constant will become zero and the equation will 
be satisfied. Thus we have scanned over different values of the boundary values till (4.10) 
evaluates to zero. The result of this is an exact solution that satisfies all the equations of 
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Figure 7: Plot of the charged W boson - numerical. The second plot expands the small u region. 
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Figure 8: Plot of (f)^ - numerical. The second plot expands the small u region. 
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Figure 9: Comparison of analytical (dashed) and numerical solutions for u < 100. The agreement is good 
for large u as expected. 
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Figure 10: Comparison of analytical (dashed) and numerical solutions for asymptotic values of u. The dashed 
line is not visible as the agreement is good. 



motion. This is shown in Figures 5,6,7 and 8. The value of the SU{2) fields for this solution 
are very small in the entire region and hence the Yang-Mills approximation to DBI is reliable. 
This then establishes unambiguously that a condensate is indeed formed. 

5. Phases as a function of fi and T 

This section summarizes the phase structure (at zero temperature) based on the calcu- 
lations of the earlier sections. To study the thermodynamics we go to Euclidean section by 
setting it = Ie where Ie has a period /3 = ^ . To compare the phases at zero temperature we 
consider the Euclidean action with /3 = oo, which effectively means consider the coefficient 
of ^. We have a background field = At. Wick rotation converts this to Atj^ which satis- 
fies Atj^dtE — Afdt. Thus At^^ — 4*- . When one does the Euclidean functional integral we 
analytically continue to real and integrate over real field configurations. However 
if there is a backgroimd real At (that is a solution of the classical Minkowski equations with 
a real charge), this has to be put in as an imaginary background value for Atj^ . Thus if 
one wants to perform the functional integral at one loop, then one has to work with real 
Atj^ fiuctuating about an imaginary background value. However at the classical level, we 
need not do all this - we can just plug in the Minkowski solution into the action without 
introducing Atj^ . 

Since the configurations depend only on u, we need only worry about the u integral. 
Thus we need to look at S[Al{u), t{u)] evaluated on the solutions. 

As a first step we just consider the solutions of the Dirac-Born-Infeld equations considered 
in Section 1 where the sources are localized as S{u — Uc)- This can tell us whether the phase 
with number being non zero is preferred over the phase with net number being zero. We 
remind the reader that we are referring to the component of the U{1) x SU (2) . [ta are the 
SU (2) generators.) We assume that the U (1) b baryon charge is already non zero and we have 
the cusp configuration. When ^3 number is non zero, one can further ask whether there is a 
^3 number violating condensate. This has to be asked within the Yang-Mills approximation. 
This is the second step. 



26 



The Dirac-Born Infeld action with source is (see equation (2.4)) 



Sdbi ^-f j da^ [u'' + u'fr'^ - /(9,ao)^]^ - UBTaoiu,) (5.1) 
and the solution for qq is 

- fiu^ + dl) ^^-^^ 

This determines aQ{u) in terms of do- In our case we have two branes. We write di,d2 
corresponding to the two branes as: 

J J , ^3 , , dz 
dl = do + y 02 = do - y 

where ^3 is the charge that we are concerned with. The chemical potential — Aq(oo). 
We will fix a gauge by setting Aq(mc) = .^^ 

We have two branes with different profiles t(m), which in turn depend on the constants 
d and c, see equation (2.10) . We need to impose that Uc is the same for both branes. So if 
di,Ci are given, then Uc is fixed and so C2 is fixed once ^2 is. Thus wc Lcgcndre transform 
r'(M) in terms of cq. Wc can also substitute the solution for aQ{u) and t{u). Thus we find 
that the Legendre transformed (with respect to r) is 



2 



nil,, Co) = ^ , ^ ° ^ , Dido, Co) = u'f + dlu'f - cl 
Dido, Co) 



Generalizing to two branes one finds: 



fu^ - cj fu^ - c 



nii^o, 1^3, Co) = du { — + 



fu2 ^yDidl,Cl) fu2 ^JD{d2,C2) 

where di = do + ^ d2 = do — ^ • 

One can compare this action with the action with da = and the same value of /X3 at the 
boundary. If da = then the branes are right on top of each other and there is no separation. 
Thus not only is there no charge, there is also no number violating condensate. We 
have an unbroken SU (2) gauge symmetry. If da = then Aq is constant everywhere and is 
pure gauge. The action is 



^o(Ato,/^3,Co) = 2 du 



fu^ - ci 



fu2 y/D{do,Co) 



^^Normally when the gauge field extends ah the way to the horizon, one usually sets the gauge field to 
zero at the horizon. The gauge field here lives on the D8 brane, which does not extend beyond u = Uc ■ 
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The precise relation between and do can be obtained by first calculating the free energy 
F{dQ, Co) with fixed dg (Legendre transform of ft with respect to // , or equivalently, of S 
with respect to Oq ). 



F(do, Co) = T / du — 3 H — Uc . 

We have included the D4-brane action of the U{1) charges. This gives 

OF ~ p , ^ f r , 'kill fu, 

A* = TTT' = T du Oudo + = T du 



ddo Juc 3 J^^ y^D{do, Co) 3 

The free energy difference is Q{di) + ^2(^2) — 2 f2((io) • This has to be plotted for different 
/X3 keeping fixed /xq = /xi + /i2- It is simpler to plot the free energy Q as a function of ^3 and 
give a separate plot giving /i^ versus ^3. 



The chemical potential is given by: 



A«3 = A*i - A«2 



^ I du { 

Juc \ 



3 3 

diU2 d2U2 



Figure 12 shows a plot of [x^ vs ^3 for the separated D8 brane solution (i.e. with a 
number) for a sample value of Cq, c?o • The free energies are plotted as a function of ^3 in 
Figure 11. Wc sec that for all values of ^3 the free energy difference is negative. For e.g. for 
d^ = 10 (so that d^ « do), we find that Afl ~ —89.97 . Thus one infers that as /X3 increases 
so does ^3 (see Figure 12), and so (see Figure 11) one has a phase transition to the phase 
with non zero number. 

We have thus established that when there is a non zero chemical potential for the 
phase with a finite number density of t-i, charge is favored. We also have exact (numerical) 
solutions for this phase which have charged fields condensing. A calculation of the free energy 
difference for the given exact solution gives ~ —4.40 * 10"^ and shows that this condensate 
is favored over the trivial solution (with condensate being zero). Thus putting these results 
together we have established that a BCS condensate is formed and we have a holographic 
description of the ground state. This is the main conclusion of this paper. We now discuss 
below the implications of this for the boundary theory. 

Our experience with the boundary BCS theory says that there should be a condensate, 
since when there is a finite number density of fermions, and hence a Fermi Surface, one 
expects the BCS instability. This should thus reficct itself in a t-i, number violating condensate 
of charged fields. As mentioned above the existence of this condensate is proved by the 
holographic calculation. 

The nature of the boundary condensate also needs to be determined. As discussed in 
the Introduction the most plausible candidate seems to be a BCS like Cooper pair BB 
if N(. is odd. We are at zero temperature and in the confining phase in the boundary. 
B = eai,a2,-,aNc''^"'^d"'^''^"'^d°'^ ' ' IS cvcu. One could have scalar baryons with 
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^3 = +1 condensate. This would be analogous to a BEC because the particles are strongly 
bound already, before they condense. 

The effect of this condensate should be visible in the conductivity. Thus if one calculates 
the current current correlator one should find a gap in the imaginary part (equivalently in 
the real part of the conductivity). This is a problem for the future. In the next section 
we take a first step in this direction by considering the conductivity of the Sakai-Sugimoto 
model without any SU (2) breaking. 



6. Conductivity 

Consider the DBI action given by equation (2.4) with ub — ^ and a — u , 

I' - 

Sdbi ^-f j d'x du^ [1 + u'fr" - fid^aof]-^ . (6.1) 
The equations of motion are 



l + u^Pr'^-fid^ao) 



2 



= Cn 



2 



dt (6.2) 



where Cq and do are constants. These equations imply, see equations (2.10), 

0,a,r = ^ , {9,rf = :^ (6.3) 
where D — u^f + d^u^f — Cq and, hence, that 

>C = ^ [1 + u'fr'^ - fid^aon ^ - ^ • (6.4) 

The D8 brane ends at u — uq where duT oo i. e. D{uq) — . This gives uq in terms 
of other parameters as 

{ul + dl)ulU = cl , h = f{uo) . (6.5) 

6.1. Fluctuaiion about the background 

Let us turn on the U{1) gauge field 02(0"^) along the a"^ = x direction of the D% brane 
worldvolume. We will consider this as small fluctuation around the background described in 
previous section and neglect any back reaction. For simplicity let us first assume 02 to be a 
function of t and u only. We will also assume a2{t, u) — a{u)e'~'^^*. Let a dot (') and a prime 
(') denote t and u derivatives. 
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For matrices M and 5M, where 5M is assumed to be small, we have 



^^^^ - ..i^m^U^M-^sM^y 



--tr{M-^6Mf + 0{SMf . 



Applying this relation to the Lagrangian Cdbi = \J det (G + F + 5F) gives 



(6.6) 



(6.7) 



where C = a/ det {G + F) . For the problem considered here, C is given explicitly in equation 
(6.4), and we have tr{{G + F)~^6F) = . Wc can construct 3x3 matrices out of (G + F) 
and 6F to compute the traces, as 6F is of the form. 



/ d2 \ 

a'2 
5F ^ -02 -^2 



I 0/ 

Let us call the truncated blocks of {G + F) and SF as B and SB respectively. 



/ 



B 



SB 



3 

-M2 



-a: 







^0 4-(i + «VV(«)2) 




\ 



3 
y 



Then, 



1 



tr{{G + F)-'5Fr = _ ((1 + KVV2)(a2)^ - u^fia',)') 



Then the Lagrangian Cf for the fluctuations is given by 

1 / {u' + dD .0 ID 



(6.8) 



(6.9) 



(6.10) 



(6.11) 



The equation of motion is then given by 




+ d^^ 



3 ^2 



a2 = 



(6.12) 
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Setting a2{t,u) — a{u)e we get 



(Ma]\^^.^^0. (6.13) 

The equation of motion is same for the world volume fluctuations considered above for 
both D8 and D8 branes. In the absence of charges, continuity of flux would require that the 
fields continue in the same direction along the other brane and reemerge on the boundary. 
So we can consider a{u) to be the field on both D8 and D8 branes, but u is a bad co- 
ordinate choice for such a representation. We can consider a change of co-ordinate given by 
y"^ — 1 — ^. y e (0, 1) corresponds to brane and y e (0, —1) anti-brane. The brane and 
anti-brane are joined at y = , the y = ±1 are corresponds to the intersection of the brane 
and anti brane with DA . We will call y = 1 boundary and y = —1 the "horizon". The 
differential equation will be solved with "in-going boundary condition at the horizon" , and 
we will use "AdS/CFT correspondence" to evaluate the boundary Green's function. 

This choice of boundary condition does not have a rigorous justification. However we 
can argue heuristically. We want boundary conditions corresponding to the normal modes 
or more correctly quasi normal modes. The eigenvalues then give the poles of the Green 
function. We motivate boundary conditions as follows. In the case of the black hole, in the 
AdS/GFT context, these boundary conditions were first introduced in [48, 49] where it was 
argued that at the horizon, ingoing boundary conditions are useful because with the usual 
Dirichlet/ Neumann boundary condition the eigenvalues are strictly real. The latter would 
give real poles for the boundary Green function and would not describe thermalization which 
requires an imaginary part for the poles. Ingoing boundary condition is one option that does 
give an imaginary part to the pole location. Furthermore ingoing boundary condition has 
the reasonable physical interpretation of complete absorption by the black hole. 

In the present case also we only have a few reasonable boundary conditions and we pick 
one that seems to give the right pole structure. Ingoing, is one such and corresponds to 
perfect absorption. Are there situations other than a black hole horizon where one expects 
complete absorption? The answer is yes. Consider electromagnetic waves in a wave guide 
with an oscillator source at one end and the other end open and radiating via an antenna 
into three dimensional space. The cross section of the wave guide is two dimensional. It is 
well known that when there is impedance matching the entire output of the wave guide is 
radiated and there is no refiection at the boundary. The radiating antenna behaves like a 
perfect absorber at the end of the wave guide. The correct boundary condition inside the 
wave guide in such a situation is the ingoing one. Thus perfect absorption is not unknown 
outside the black hole context. In our case we have a lower dimensional wave propagation 
channel opening out into a higher dimensional space - the cross section of the D8 brane is 
effectively three dimensional and is connecting to a four dimensional D4 brane. This gives 
a possible physical motivation for considering ingoing boundary conditions in the present 
situation. 

We will set ukk = 1 in all the future calculations. The differential equation for a{y) near 
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"horizon", i.e. near y — —1 , is given by 



, (Pa 1 da 2cu^ ^ ,x 

which gives the near-horizon behavior of a{y) as, 

a ~ e ^ (6.15) 

The solution a e ^ v^T+y gjygg ^j^g in-going wave. Let a,{y) = e * ^ ^^^^p{y) ■ Then the 
differential equation for p{y) is solved numerically to obtain a{y), with boundary condition 
that p{yh) = 1 + ^^Vl + Vh and p'(i//,) = , where yn is a cut-off near y = -1 . 

The expression for p{yh) is obtained by doing a series solution of the differential equation 
near y = —1 . 

The differential equation for a{y) near "boundary", i.e. near y = , is given by 

, , (fa 1 da 2u'^ ^ . . 

The solution for a{y) near boundary, i.e. near y — +1 , behaves as 

a{y) ~ A{l + A^{l-y) + A2{l-yf + ---) 

+B (1 - y)l (1 + - + ^2(1 -!/)' + •••) (6.17) 

where A and i? are arbitrary constants determined by boundary conditions. Ai and Bi are 
constants determined in terms of parameters of the differential equation. So, the Green's 
function is given by, 

GR{uj,Uo,do) ~ — = hni ^ — . 6.18 

A y^i a{y) dy^ 

The conductivity can then be obtained by 

(7 = — . (6.19) 

Figures 13 and 14 show variation of conductivity with frequency. We have set do = in 
our conductivity calculations as the variation of conductivity with do is very small. 

6.2. Modified Sakai-Sugimoto model 

In the original paper, the scalar bi-fundamental tachyon field arising from open string 
between DS and D8-brane was neglected. Later various authors [27, 34] have included 
the contribution of the tachyon field which is expected to be the origin of quark mass and 
condensate. The effect of tachyon field for the case of gauge field fiuctuation is that the 
fiuctuations become massive, with mass proportional to the tachyon field. The tachyon 
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Figure 13: Real part of Conductivity with do = 0: uq = 10 (Black, continuous) and uq = 15 (Red, dashed) 



Figure 14: Imaginary part of Conductivity with c?o = 0: uq = 10 (Black, continuous) and uq = 15 (Red, 
dashed). 



field has an approximate behavior like e "o asymptotically. We will assume the mass term 
corresponding to the a2 fluctuation is proportional to mle "o , to get a qualitative nature 
of the effect of tachyon on the conductivity. More rigorous analysis needs to be done for 
quantitative results. Let us consider the modiflcation of equation (6.13) given by 



A similar analysis, as in the previous section, can be done to obtain conductivity. Figures 
15 and 16 show the conductivity for various values of mo . We see that with a non-zero 
value of mo the conductivity develops a mass-gap. Effect of do on the conductivity is very 
small, and we have set do = for all our numerics. As we increase mo, both mass gap and 
peak heights increases. Also, for any flnite value of mo, a pole appears at w = for the 
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-10 





(6.20) 
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imaginary part of conductivity (fig. 16), which imphes a delta function peak for real part of 
the conductivity by Kramers-Kronig relation. The delta function function peak can not be 
captured by numerics. 

Re(o-) 

; : 150 ^ : ; 




Figure 15: Real part of Conductivity with do = 0, uq = 10: toq = (Black, continuous), mg = 5 (Red, 
dashed) and toq = 10 (Blue, dotted). As we increase the mass mo, the peak height increases. The peak for 
mo = 10 is outside the range of the ordinate value shown in the figure. 




Figure 16: Imaginary part of Conductivity with do = 0, uq = 10: mo = (Black, continuous), mo — 5 (Red, 
dashed) and mo = 10 (Blue, dotted). The peak for mp = 10 is outside the range of the ordinate value shown 
in the figure 



7. Summary and Conclusions 

In this paper we have attempted to find a holographic description of the BCS Cooper 
pairing instability and thence describe the strong coupling microscopic BCS phenomenon. 
The analogy with chiral symmetry breaking in QCD was exploited and the Sakai Sugimoto 
model with finite number density of flavored fermions was studied. The background charge 
corresponds to a U{1)3 embedded in the SU{2) flavour- The flavor SU{2) was broken to U{1)3 
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by the background and it was shown that in this situation there is a tachyonic instabihty 
that causes fields charged under this f/(l)3 to condense. This is the bulk description of 
the BCS instability. Analytic solutions describing this condensate were given in flat space. 
The stability properties were also studied. The curved space equations were then solved 
analytically in some approximation. The exact solutions have also been done numerically. 
The free energy difference shows that this solution is favored over the trivial solution. This 
solution then describes the BCS ground state of the superconductor. 

Once we have this solution one should be able to calculate various quantities such as the 
AC conductivity using the AdS/CFT dictionary. However there are two complications. One 
is that the bulk is not asymptotically AdS. The second is that the flavor branes do not extend 
all the way to the interior. Thus the usual "infalling" boundary conditions prescription has 
to be modifled. A preliminary exploration of this problem in the original Sakai-Sugimoto 
model was done in Section 7. The chiral current Greens function was calculated. It has the 
expected properties of a mass gap. 

The detailed analysis of this in the context of the Sakai-Sugimoto model as well as in the 
modifled form in this paper need to be studied. Another important question is to look for 
the bulk signature of the Fermi surface in the boundary. 

We hope to return to these questions soon. 

Acknowledgements: We thank G. Baskaran and R. Shankar for many useful discus- 
sions. We also thank the unknown referee of Nuclear Physics B for very useful comments. 
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Appendix A. D8 brane with non zero electric field 

In this Appendix, we formulate the equations describing the profile of D8 brane in the 
background of D4 branes, and with a non zero world volume electric field turned on. 

Wc write the Dirac-Born-Infeld part of the D8 brane action in the ten dimensional 
background of Di branes as 

Sdbi = -Cg J d^cr 9se~^ V -det G 
where the components G^i, of the matrix G are given by 

Giiv — 9mn d^X^ di,X^ + F^i, 
with qmn being the ten dimensional background metric, 

(a;^«,a;^a;^a;^r,^^■■■ ,^^) 

are the ten dimensional coordinates with r ~ r + 27rr4, and F^^^ = d^ai, — d^a^ is the U{1) 
gauge field strength on the brane. We choose the nine dimensional brane coordinates cr'' to 
be given by 

Wc consider the case where is the only non vanishing component of the U{1) gauge 
field and (m, r, a^) depend on a only. The functions u{a) and T(cr) describe the DS brane 
profile in the [u, r) plane. It then follows that 

Gaa — 9uu IJ-a + 9tt ' ^Oa — ~GcQ — Foa — ~{0'o)a , 

and Gfj,i, = g^i, otherwise. The subscripts a on {u, r, ao) denote their 
cr— derivatives. 

In order to solve the equations of motion which follow from the action Sdbi and to obtain 
the cr— dependence of {u, r, ao), note that Sdbi now resembles a world line action 



'DB7 



do- \l i^ia 



where = (m,t, Oq), the subscript cr on ^ denotes its cr— derivative. For our case, qmn is 
diagonal and we have = diag {hi, /i2, h^) where 

hi = gle'"^^ (- det g^^) \a=u 
hs = 

and depend on = u only. The consequent 'geodesic' equations then give 

h2 r„ = Co (A.l) 

hs {ao)a = do (A.2) 

hi ul + h2 tI + {ao)l = E (A.3) 

2 hi + {hi)u ul - {h2)u rl - {h^)u {ao)l - (A.4) 
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where Cq, do, and E are constants and the subscripts u on the hs denote their derivatives. 
It follows from the above equations that 



The right hand side of the above equation is a function of u . Generically the location 
u — uq oi its zero, if exists, denotes a turning point for u{a). Thus if u{a) — )■ uq from above 
as o" — > (To then, generically, u starts increasing as a increases beyond uo which can be seen 
from Ufja equation. Also, generically, the evolution of t{(j) and 00(0") is monotonous across 
(To . In the context of D8 branes in the background of DA branes, uq will denote the lowest 
point of D8 brane profile in the (m,t) plane. Also, {ao)a equation above with do ^ shows 
that such a D8 brane can support a non trivial U{1) gauge field on its world volume. It is 
now straightforward to obtain the evolution in terms of u alone using = ^ and similarly 
for (ao)u . Using Fq^ — — (ao)w a-nd the expressions for /i2 and in equation (A. 3), we 
have 

A ^ 1 + g^-g^r tI + g^'g^- F^ = (A.6) 

where is given in equation (A.5). 

Writing the ten dimensional background fields of DA branes as 

d'SiQ = gMNdx^^ dx'^ 

(A.7) 

and e* = gs where H = ^ , and f — 1 — , we have 



where D = {Eu^ + §)f - cl , and 



j2 „2 TT 

2 «o 1 ^0 ^ 



l«o;„ - ^ ^ ' - /2 L> • 

The expressions for (ao)^ and given in equations (2.10) correspond to the choice E — \ . 



Appendix B. Embedding of the SIJ{2) group in SU{2Nd) 

We have two sets of coincident branes. For concreteness and clarity of presentation let 
us assume that each coincident set has three branes, i.e A''^ = 3. (The minimum required 
is two). So we have six branes in all. Our configuration breaks U{1) x SU{6) down to 
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U{1) X U{1) X SU{3) X SU{3). We have 36 generators in U{6). We can write them using 
the direct product notation. Let A stand for the Gell-Mann Lambda matrices for the SU (3) 
part: A^'^' -^ and A° be the identity matrix. 



Similarly a^'^'^ are the SU{2) matrices with being the identity. 

Then the unbroken generators of U (6) arc: A'^ (g) a^, A° (g) are the two unbroken U (l)'s. 
Schematically (the labels 1,2 stand for the two branes): 

or equivalently: 



y ' - \ AO 



V{\)b^{^^1 °o); t/(l)3 : 2^3 = ( _°o) (B.l) 

Similarly A" ® {a^ + a^)l2, A" ® ((7° - a^)l2, a = 1, 2. ..8 are the two SV{3). These are 
not excited in our condensate solutions. They do contribute in an instanton contribution 
localized u — Uc- They generate the Chern Simons interaction which gives a source at 
u — Uc- Schematically: 

■ A" A / 



; ' V --^^ ' ^^'^^ 

The eighteen broken generators are: 

A° (8) (T^'^ are the generators corresponding to the VF^'^ massive gauge fields of the SU (2) 
that condense. Schematically 

^1 = ( °o ) ; ^2 = ( iAO 0^ ) ^^-^^ 

Finally the sixteen A" ® a^'^ are the broken generators corresponding to gauge fields that 
are not excited. Schematically: 



A" y (0 -iX" 
A" J ' I iA" 



(B.4) 



The total number of generators is 36 as required for U{6). This is easily generalized to 

U{2Nd). 

The Yang- Mills analysis in this paper involves ti,t2,^3- As can be seen they see nothing 
of the SU (3) internal structure. Hence these indices are suppressed in all the equations of 
this paper. 

Thus for each set of three coincident D8 branes the DEI analysis given in this paper 
follows that of [20, 34]. 
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Appendix C. Symmetrized Trace Prescription 



In this appendix we give a useful result on symmetrized trace (Str) prescription. This is 
useful in the two flavor case when there is only the component of the U (2) Yang-Mills fields 
being non-zero and the branes are on top of each other (so that the symmetric "metric" part 
of the Dirac-Born-Infeld is proportional to the identity). Evaluate using symmetrized trace 
prescription the object 



The prescription is to expand in a power series in F3 -|- F and symmetrize each term. 
Let F = F+T^ 
Observations: 



Let F = F+T+ + F'T- and F3 = i F\^ 



1. Since F^o = — -fo„ is the only non zero component, the trace (over Lorentz indices) of 
an odd number of F's is zero. 

2. 



^3 4 

p2 ^ p+p-^^+^^-^ ^ F+F'I 
F3F + FF3 = . 

3. Consider a term F^F^. It has to be symmetrized completely so we have all per- 
mutations. Consider the first two places: If it is F3F..... then one has to add the 
permutation FF^.... where the three dots are identical in both cases. This will give 
F3F... + FF^... = {F3, F}... = . Thus we can conclude that in the first position either 
there should be two F's or two F3's. Having done this operation on the first two places, 
we consider the next two places. The same argument holds. This can be repeated. 
Thus we conclude that Str {F-^F"") ^ (F2)f (F|)t . (Here F| = ^ and F^ = 
F~^F~, i.e. without any matrices.) We need to determine the precise coefficient. 
There are ^"t",-*' permutations. So 

T). 777. ^ 



Str (F^F"") = ———7 Tr V (Fg^F™) . 
^ ^ (n + m)\ ^ 

Perm 

The number of non zero terms is the number of ways we can pick y places to place 



m-{-n I 



2 2 



the F' in a total of ™ places (i.e of the form F^F^F^F'F^... ). This is nfer- So 
the final answer is (A factor of 2 for the trace) 

m+n I 
I 2 • 

'f!f! (n + m)! 
Based on the above n, m have to be even. 
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Our starting point is 



Exp 



^ 2 + 4 + 6 + 



using the fact that odd powers vanish. This is the power series expansion of ^/l — (F + F^y . 
So let us define C„ by: 

oo 
n=0 

Then 

n 

Thus the coefficient of a term FnF^ in the above is Cn+m I^OllH^ _ Thus using the formula 

•1 — 2 — n\m\ ° 

for symmetrized trace derived above we get: 

(n + m)! ^! n!m! 2^^/p2^g 

m+n I 



2 ■ 2 ■ 

Wc can now use it to get coefficient of powers of : 
Case 1 m = : This gives 

2 (F|)t . 



Summing over n (even) gives 2yT^-^F^ 



Case 2 m = 2 : 



2C^t+i (^ + 1) 



We need to perform the sum (write n = 2ni) 

oo 

2F^5]C„,+i(ni + l)(F32y 



ni=0 

Let us rescale F^ xF^ . Then wc get: 

2F2^C„,+i(ni + lK^(F32r 
m 

which can be rewritten as 



dx 

m 



— 9 ,^ni+lrF2\ni+l 
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n2=l 



where, in the last expression, n2 — rii + 1 . We can extend the lower limit of sum to 712 = 
by the addition of an x-independent term without any change in it's value. Thus we get 



n2=0 



F 



2 



(setting X — 1). 

Thus the first two terms in the power series for is 



2^ 1-F. 



'2 



General Case: 

Let us consider the general term with Let | = ni and y = mi. Then we have 



2a 



ni!mi! 



o/p2vmi^ K + + mi - + 1) 

= ^l-f" j W+mi j [^3 

mil 

Replacing — > xF^ as before we can write this as: 

1 

9/7Ti2\mi^ ™ni+mi / 7Ti2\ni 

2(F) ^n,+m,:^J^X (F3) 

p2 1 ^mi 
— Of ^m+nii/ 77i2\ni+mi 



Do the sum over ni. 



.F2_ 1 d'"! 



711=0 



Let ni + rrii — n2- Then we get: 



p2 1 jmi °° 

^Ff mi!(ia;™i ^ ^ 



n2=mi 



42 



Because of the derivative we can extend the sum over ^2 to 0. 



oo 



3 ^ 712=0 



Summing over mi one finds a Taylor series for 



l-a;F| . (C.2) 



2^^l-^-^-F+F- . (C.3) 

Actually one can conclude that this had to be so from the first term in the Taylor series, 
(the case m = 2) and requiring Lorentz invariance of the final expression, which uniquely 
fixes the expression inside the square root. 

Appendix D. Nf number of DS branes: Action 

In this Appendix, we consider the non abelian generalization of the action Sdbi for Nf 
number of D8 branes. We set u = it . Thus, the worldvolume coordinates cr^ are 

The D8 brane action may now be written as 

Sdbi = -Cg J d^a QsC'^ Str -det {G + M) 

where the components G^i, and M^u of the matrices G and M are given by 

G^u = Qtiu + Qtt dij,T dvT + F^j, 

with qmn the background metric of DA branes. In the above expressions, we have spUt the 
U{Nf) terms into the abehan U{1) terms and F^j, , and into the non abehan SU{Nf) 
terms V^cf) — {T>^(f)"-) ta , (f) — cf)"' ta , and f^i, — /^^ ta where ta are the SU{Nf) generators 
satisfying the algebra [tb, tc] — i ta and 

= d^Al-d,Al+ gC\,AlAl . (D.l) 

The SU (Nf) invariant action is to be obtained by symmetrized trace (Str) prescription. In 
this paper, g — —1 . 
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In this paper, we will keep r and F^i, terms fully, but keep the non abelian and f^iy 
terms only up to quadratic order. We have 

_ 1 Q^u^xa M,xM^^ + 0(M=^) (D.2) 
where G^^ are the components of the matrix . Let 

denote the symmetric and antisymmetric part of G^^ . Taking the symmetrized trace over 
the group indices we then have 

G^'^M,^ = G^^ grrV^rT^A"" 

up to quadratic order in 0" and /^j, . The negative signs in the last expression arise because 
of the shuffling of i/, • • • indices. We write the resulting action as 



Sdbi — — Gg 



Jcfa-fC^ , 7 = t,,e-* V-det G (D.3) 



with given, up to quadratic order in (p"" and , by 



+ \ (Gf G^'^ - GrG^'^) . (D.4) 

The equations of motion for the fields (0", A^) can now be obtained. The G_ terms 
cancel each other if the rank of the matrix G_ is two (which is the present case) or three. 
Ignoring therefore the G_ terms, the equations of motion are given by 

V, (7 Gf g^r ) = (D.5) 
V, (7 Gf Gi<^ = g GV 0" (7 g^r V,ct>'^) (D.6) 

where V,{*^) = d,{*^) + g C\ (*^) . 

In our case, oq is the only non vanishing component of the U{1) gauge field and (r, oq) 
depend on u only. Hence Fq^ — — (ao)u , 

Guu = Quu + 9tt Tu ) ^Ou — —GuO = Fqu , 
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and = g^i, otherwise. The subscripts u on (r, Oq) denote their 
li— derivatives. G^'^ and G'^ are now given by 

„00 „uw 

„00 „uu 

and G^'^ = G^ ~ g^^ and G^ — otherwise. A in these equations is defined in equation 
(A. 6), namely 

A = 1 + g^'^g,, tI + . 
With these substitutions, and noting that 

7 = ^,e-* \/-det G = ^ \fK (D.8) 

for the DA brane background given in equation (A. 7), it can be checked that the action Sdbi 
above gives the action in equation (4.4). For the DA brane background, we also have 



7 = 



D f 



where D — {u^ + — cl and the terms inside the brackets above — )■ 1 for large u . 

The analysis leading to the action given in (D.3) is also apphcable to the case where, 
instead of non abelian fields, one switches on other component(s) of the abclian U{1) field 
and studies their leading order fluctuations. For example, let a2{t,u) be the U{1) gauge 
field along the = direction in the D8 brane world volume. Then the action is given, 
up to quadratic order in 02 , by 

S{a2) - y" 7 G^^Gl^ {d,a^ - d^a^) {d^a, - d,a,) 

~ y 7 (Gf G^2 {dta^f + G'T G^2 {dua^f) . (D.IO) 

The second line follows since 02 is a function of (cr°,(T^) = {t,u) only. Using equation (D.9) 
and G'^ — g'^'^ — y/H , we get the action given in (6.11) 

S{a2)^ I VhD (^-^^^^^ {dta^r + {d^a^r^ 
since ^^^^ - {u' + dl) . 
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Appendix E. Nf (= 2) number of D8 branes: 
Equations of motion 

In this Appendix, we write down the general equations of motion which follow from the 
action given in equation (D.3). These equations may then be specialized to the various cases 
studied in the paper. 

Consider the case where ^4^ = if /U ^ 0, 1 and the fields {Aq, A^, 0") depend only on 

----- - - - - - - - 



((jO, a^) = {t, u) . Then /^^ = if {fx, u} 7^ {0, 1} . Taking G^^ to be diagonal and defining 

15 

7"^ = 7 grr , 7' = 7 Qrr , 7^ = 7 G'^GT 
the equations of motion (D.5) and (D.6) for the fields (0", Ag, A") become 

Po(7°W) + 1^1(7' ^^10") = (E.l) 
1^oi¥ /oi) - g 0' (7' ^^10^) = (E.2) 
(7^/01) + 5 / (7° ^^00^=) = . (E.3) 

We now specialize to the case of SU{2) , namely Nf = 2 . We choose the SU{2) generators 
to be t_, ts) with the independent non vanishing components of the structure constants 
C\c = -C''eb given by 

G\_^_ — — i , G ^_ — + i , C^+_ — — 2i . 
The fields (0", ^4^) are given by 

= (6^^ t„ = 0+ t+ + 0- t_ + 0^ ta 

where (0'^, are complex fields, {(t)~,A~) their complex conjugates, and (0^, A^) are 
real fields. Then V^j^cf)"- are given by 

P^0± = T ^9 (40^ - 40') (E.4) 



and/^^ by 



= - 2^(7(A+r-^^0+) (E.5) 



f% = d,A^ - T - A^^Al) (E.6) 

41 = ^M^' - d.Al - 2ig {A+AZ - A' At) . (E.7) 



^^For the DA brane background, using grr = and equations (D.9), we have 
Note that 7^ » |7°| ~ |7-^| for large u since H = . 
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We further choose a gauge where Aq = = . (A" corresponds to the gauge field along 
= u direction.) The non vanishing gauge field components are then [A^, A"^) = [A, w^) . 
Then 



(Po0^, Po0') = {(t>fTigA(t>^, 0?) (E.8) 

(E.9) 

ifoi, /oi) = -AJ (E.IO) 

where the subscripts t,u on the fields w^) denote their derivatives with respect to 

t, u . 

The equations of motion (D.5) and (D.6) for a = 3 become 

(7V?)* + {l'4>l)u - 2tg{^\w^<p- - w-<P^))u 

-2ig'y\w+(j)- -w-(j)+) -4g^'y^w+w-(j)^ = (E.ll) 

+25V(w+0" = (E.12) 

-(7^A)« - 2ig-ff{w+w^ - W-W+) - 2zgj%(f)+(f); - (/)-(/)+) 

+V(7-^w+w- +7V^0")^ = , (E.13) 

and, for a = + , they become 

(7°(0i+ - tgA<p+)), + (7^(0+ + tgw^<p'))u - tg^'AcPt 
+ig-f^w+(pl + 2/7^(w+0~ - w~(j)+)w+ - g'^-f'^A^(j)+ = 



(7^(«;+ - igAw+))t - igi^ Aw^ + igi\(t>^(t>t - (t>^(t>l) 
-2g''^\w+(t>- - w-(t>+)(t>-^ - /(7'(0')' + ^^A^)w+ = 



(E.14) 



(E.15) 



{l^wt)u - igi\(t>^(t>t - <P^<f>t) - ig{l^w+A)u 

-ig-ffw+Au-gyA(l)^(l)+ = , (E.16) 

and the a = — equations of motion are the complex conjugates of the a = + ones. In the 
static case, all derivatives vanish. Let (f)^ — (j) e^^^ and — w e^'^^ where (f) and w are 
real. It then follows from equations (E.12) and (E.15) that Q, — 9 -\- — constant . 

In this paper, g = —1 . Setting 7° = 7/ = — 1 and 7^ = 1 and charged scalars 
to zero, gives the fiat space equations of motion (3.3) - (3.7). With 7*^ = yj—gg^^grr , 
7^ = \/—gg^^gTT , and 7^^ = ^/—gg'^^g^^ , and with all t— derivatives vanishing, the above 
equations give static curved space equations of motions (4.5). 
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Appendix F. Flat Space Equations 

Two solutions were given in the text. The equations of motion are given here. The 
field X was defined in Section 3. If below we assume that x = 0i is real then we have the 
solution with one adjoint scalar corresponding to the D8 brane case. If we leave x complex 

the equations could be used for other branes. 

The equations of motion, in the same gauge as above, and with time derivatives set to 
zero, are: 

-dlAl + 2 {2w^w- + x+X" + X-X^) Al (F.l) 



- x^duX - X duX^ + X'^duX - X duX^ 



+iw+ (x'x~ + X'X") + (x'x"" + X'X+) (F-2) 

(5S 1 

T~3 = n^uX^ + 2^ (^^^«X" - w'duX^) + i ix'duW^ - X^duW) 
OX 2 

-\ [x'itx' - X'X ) + X"(x'x+ - x'x+)] - 2«;+«;- t (F.3) 
: 2i d^{w^Al) + 2i w^d^Al - ix'x'- + tx") Al (F.4) 



5Ao 



= i ix-'d.x' - x'duX-") + i ix'-duf - fd^X^) 



+ 2 {{Alf - x+X~ - X'X"" - XY) w+ + 4 X+X+ (F.5) 



(5S 

~- 9lx^ - 2i w'^duX^ - i x^duW^ + 2 (w+x~ - w~x^) w+ 



^X 

+ {Alf x+ + ^ x'(xV - X'X+) - X+(X+X~ - X"X+) (F.6) 

Of course if we set x = wc have the same solution as before. Furthermore, if we set 
one of the scalar fields, such as 0i = , then the system is again the earlier one and we only 
have an analj^ic solution where 02 is also zero. So we try to set a different set to zero. One 
can try for instance to set (f)f = = (f)f . 

The Aq equation in fact suggests this possibility: x'^ = X"^ x^ + X^ = • Thus the 
term (x^X''' + X^X''')^o vanishes and we get the same relation as earlier between A^ and . 
Note that x+ = (x")* • Thus x+ = -(x")* ■ 

In the Al^ equation, the coefficient of w vanishes and the derivative terms become: 
2x~^«X^ — 2x'^duX~ ■ This vanishes if the phase of x~^ is constant, independent of u . 

The equation reduces to (Aq)^ = (x^)^ if we choose x^'W~ = X~'W^ 5 which means 
X~^w~ should be imaginary. Thus we let X^ = '^X = ~X~ '^^^ = = w , so the phase 
of w'^ is also It- independent. 
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Using x~^'^~ — X~'^~^ ^-iid the fact that the phases of x ^-nd w are constants, we see that 
in the x^ equation the derivative terms canceL And using the reahty properties of x , we 
get the same equation as ear her: 

Note that x'^X~ is negative definite. 

Finally the x~ equation reduces to dlx~^ = and 2iw~^duX^ + ix^d^W^ — . We impose 
these two separately, because has to satisfy the same equation as . 

Thus, when all the dust settles we have a system of equations, very similar to the earlier 
one, except that x"*" can be nonzero: either a constant or linear in u. But this changes 
fairly dramatically the behavior of Aq - it becomes exponential rather than linear. If we take 
X"*" = Xo = constant then A^ can be solved for in closed form: One finds using the same 
methods as earlier 

^o(m) = yaXo cosh (xo u) (F.7) 
Note that duAl{0) = , and thus the electric field vanishes at ti = as required by symmetry. 



Appendix G. Nf (=2) number of D8 branes: 
Set up for the study of Stability 

In this Appendix, we write down the general equations of motion for small field fiuctua- 
tions in the background of non zero A^ and . They may then be speciahzed to the case 
studied in the paper. 

Consider solutions to the equations (E.ll) - (E.16) for the static case and with 0=*= = 
— . Setting cf)^ — B , equations (E.ll) and (E.13) give 

{'y'Bu)u = {'y^Au)u = . (G.i) 

Hence = q and 'j'^Au = E where q and E are constants. Other equations are satisfied 
identically. Setting B{uc) = A{uc) = , we have 

Consider the fluctuations of the flelds around this static background. Thus, we write 

A{t,u) = A{u) + a{t,u) , B{t,u) = B{u) + b{t,u) 

where A{u) and B{u) are the static background solutions given above. The flelds (a, w^, 6, 0^) 
are functions of (t, u) and are assumed to be small. Their equations of motion follow from 
equations (E.ll) - (E.16). Noting that (7°, 7^, 7^^) depend only on u and that A{u) and B{u) 
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are static background solutions, we have, to the hnear order in (a,t(;^, b, 0^), 

il\)t + {lX)u = (7^««)t = (7^««)« = (G.2) 

+ig^^B^w+ - = (G.3) 

-ffw+ - 2igjfAw+ + igi\B(t)l - B^(f)+) 

-/(t'^' + 7^^')^+ = (G.4) 

-ig-ffA^w+ - g^j^AB(f)+ = . (G.5) 

In the above equations, A and B are the static background solutions and the t, u subscripts 
denote t and u derivatives. 

The independent set of fluctuations are (a) , (6) , and [w^^cj)^) . Wc have that (7'^a„) 
must be constant and, hence, it simply shifts the background constant E . The h fields are 
traveling wave type fluctuations. In the following we set a = 6 = and take (w^, 0^) to 
be given by 

^+{t, u) = 0(M)e-^"** , u) = iw{u) 

where m is a constant and {w, (p) arc functions of u only. After a little algebra, their linearized 
equations of motion given above may be written as 

(7V«)«-5((7'w)„S + 2y^/;S„)-7°(m + ^A)2 = (G.6) 

g'j' {B(l)^-B^(P)-{'yf{m + gAf + gYB^)w = (G.7) 

{{yw)^{m + gA) + 2g-ffwA^)-gj\m + gA)B4> = 0. (G.8) 
An equation for alone can now be obtained. For this purpose, let 

hA^-ff{m + gAf , he^g^B^ , h^hA + he ■ 
Then equation (G.7) and equation (G.8), multiplied by (m + gA) , may be written as 

gj {B4>^ - B^(j)) -hw = = {w hA)u r ^^^0 ■ 

T 

Using the above equations or, equivalently, equations (G.6) and (G.7), it can be shown after 
a straightforward algebra that 

7°/i , B^ / Ha 



+ ( In 



= . (G.9) 
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It can be checked that equation (G.6) follows upon using equations (G.l), (G.7), and (G.8). 
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Consider the case where g — —1 , — — —1 , and 7^ = 1 . Then B — qu , A — —Eu , 
and 

Ha = —{m + Eu)'^ , hs = (fu^ , h = (fv? — {m + Eu^ . 
Equations (G.6), (G.7), and (G.8) become 

+ {wy,B + 2wB^) + (m - A)2 = (G.IO) 
- + (S^ - (m - w = (G.ll) 

(-^^(m - A) + 2wA„) - (m - = (G.12) 

which are same as equations (3.11) - (3.12) (Amongst these three equations only two are 
independent). Equation (G.9) becomes 

2mq^u f 2mg2 \ 



(m + £Jm) /i " \ (m + £^w) /i 
since, now, 7^ = 1 and 



hJu h \ Hb)^ {m + Eu)h 

This is same as equation (3.13) for with obvious identifications. 
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